The sun's position in the sky 
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We express the position of the sun in the sky as a function of time and the observer's geographic 
coordinates. Our method is based on applying rotation matrices to vectors describing points on the 
celestial sphere. We also derive direct expressions, as functions of date of the year and geographic 
latitude, for the duration of daylight, the maximum and minimum altitudes of the sun, and the 
cardinal directions of sunrise and sunset. We discuss how to account for the eccentricity of the 
earth's orbit, the precessions of the equinoxes and the perihelion, the size of the solar disk, and 
atmospheric refraction. We illustrate these results by computing the dates of "Manhattanhenge" 
(when sunset aligns with the east-west streets on the man traffic grid for Manhattan, in New 
York City), by plotting the altitude of the sun over representative cities as a function of time, 
and by showing plots ( "analemmas" ) for the position of the sun in the sky at a given hour of the day. 

Keywords: celestial sphere, rotation matrices, calendar, equation of the center, equation of time, 
bJ[) precession, Manhattanhenge 
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I. INTRODUCTION 

Some say he bid his Angels turne ascansc 
The Poles of Earth twice ten degrees and more 
From the Suns Axle; they with labour push'd 
Oblique the Centric Globe: Som say the Sun 
Was bid turn Rcines from th' Equinoctial Rode 
Like distant breadth to Taurus with the Seav'n 
Atlantick Sisters, and the Spartan Twins 
Up to the Tropic Crab; thence down amainc 
By Leo and the Virgin and the Scales, 
As deep as Capricorne, to bring in change 
Of Seasons to each Clime; else had the Spring 
Perpetual smil'd on Earth with vcrnant Flours, 
Equal in Days and Nights, except to those 
Beyond the Polar Circles; to them Day 
Had unbenightcd shon, while the low Sun 
To recompcncc his distance, in thir sight 
Had rounded still th' Horizon, and not known 
Or East or West, which had forbid the Snow 
From cold Estotiland, and South as farr 
Beneath Magellan. 

— John Milton, Paradise Lost, book X fTl 

Growing up in a tropical country, my first interest in the problem of finding the position of the sun in the sky grew 
out of trying to determine the dates of the year when the noonday sun reaches the zenith. Later, when I first travelled 
beyond the tropics, I was fascinated by the fact (made particularly remarkable by its novelty) that in summer the 
daylight could last well into the evening. During my time in secondary school I made some attempts to characterize 
the sun's position in the sky, but lacked the mathematical tools to solve the problem in generality. 

Here I present what I have found, at length, to be the most transparent way of treating this problem rigorously. 
It is based on describing the position of the sun on the celestial sphere (a concept that should be very familiar to 
amateur astronomers) and performing several coordinate rotations on that sphere. The idea is to begin with the 
ecliptic reference frame, in which the position of the sun during the year is most easily and directly expressed, and to 
end with a terrestrial reference frame, defined with respect to an observer standing at a given point on the surface of 
the earth, at a given time. 

The mathematical training needed to understand this derivation is that which a college student should have after 
a first course in linear algebra, since coordinate rotations will be described by matrices acting on three-dimensional 
vectors. Familiarity with the transformation between rectangular ("Cartesian") and spherical coordinates will be 
helpful, but shall not be assumed. 

This work will allow us to arrive at mathematical expressions for the sun's altitude above the horizon and for its 
geographic azimuth (i.e., its compass bearing), as functions of time and geographic location. With an additional bit 
of geometry, we also obtain direct expressions, as functions of latitude and date, for the maximum and minimum 
solar altitudes, the number of continuous hours of daylight, and the cardinal directions of sunrise and sunset. We will 
illustrate these formulas by plotting them for representative cities. 

This pedagogical discussion also provides an opportunity to mention several interesting issues in celestial mechanics, 
such as Kepler's "equation of the center," and the precessions of the equinoxes and the perihelion. Another issue 
of astronomical interest that will be discussed is how the refraction of light, as it passes obliquely through the 
earth's atmosphere, affects the apparent altitude of a celestial object. We will also describe the phenomenon of 
"Manhanttanhenge," when pedestrians in the borough of Manhattan, in New York City, may see the sunset in 
between the skyscrapers. Finally, we illustrate the concepts of the "equation of time" and the analemma. 

Computer codes are readily available on the Internet to find the position of the sun in the sky (see, e.g., 0). The 
standard reference used in designing programs that compute the positions of celestial objects as functions of time 
( "ephemerides" ) is 3]. Almanacs such as [4] also provide accurate values and formulas. However, to the best of my 
knowledge there is no self-contained, analytic treatment of this problem, suitable for students without specialized 
training in astronomy or geodesy. 

Astronomical and geodetic jargon will be avoided, or confined to footnotes, except insofar as it significantly con- 
tributes to the argument's precision and clarity. Angles will generally be expressed in radians and written as dimcn- 
sionless numbers. Geographic latitudes and longitudes, as well as solar altitudes and azimuths, will also be expressed 
in degrees (identified by a superscript o) when convenient. In some cases, arc minutes (defined as sixtieths of a degree 
and identified by the symbol ') will also used. In terms of notation, the guiding concern will be to achieve as much 



FIG. 1: In spherical coordinates, a three-dimensional vector r is expressed in terms of a radial distance r, a polar angle 9, 
and an azimuthal angle <f>. The vector p is the projection of r onto the x-y plane. In terms of the rectangular coordinates, 
x = r sin 6 cos <j>, y — r sin 9 sin <f>, and z — r cos 0. 

simplicity as possible without departing too far from the established usage. 1 

The Mathematica notebook used to compute the solar altitudes and azimuths, and to produce the corresponding 
plots, is included with this arXiv submission (filename notebook/SunPosition.nb). Interested readers are encouraged 
to use this notebook to explore the derivations in this article, extending or modifying the computations as they might 
see fit. 



II. POSITION OF THE SUN IN SPHERICAL COORDINATES 

A point in three-dimensional space may be characterized by the spherical coordinates (r, 0, 4>), where r is the radial 
distance, 9 is the polar angle, and <f> is the azimuthal angle. In terms of the rectangular coordinates (x,y, z), we have 

(x \ I r sin 9 cos \ 
y = rsin9 sin</> , (1) 
z I y rcos9 J 

as illustrated in Fig. [T] 

The celestial sphere is an imaginary spherical surface, sharing a center with the earth's globe, and with a very 
large, indefinite radius. The positions of the stars, planets, and other heavenly bodies are characterized by their radial 
projection onto this surface. The largeness of the radius of the celestial sphere, compared to the radius of the earth, 
allows us, when convenient, to picture it as centered at the position of an observer standing on the earth's surface, 
rather than at the center of the earth. 

For simplicity, we take the radius r of the celestial sphere to be equal to 1 (in undetermined units) . We shall use 
a subscript (the astronomical symbol for the sun) to indicate that a vector or a coordinate thereof refers to the 
position of the sun. 



1 For instance, some simplification could be achieved by working with the geographic co-latitude (i.e., the complement of the latitude), 
but I prefer to avoid this in deference to the widespread and long-established usage. 
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FIG. 2: The sun moves along the ecliptic during the course of the year. In the ecliptic frame of reference, the sun's polar angle 
is fixed, #0 = tt/2, while the azimuthal angle 4>q increases with time at an approximately constant rate of 2n per year. 



A. Ecliptic frame 



From the earth, the sun appears to move, against the background of the distant stars, along a great circle on the 
celestial sphere called the ecliptic. 2 We will therefore start by working in an "ecliptic frame," in which the position of 
the distant stars is fixed, 3 and in which the polar angle of the sun is always #© = 7r/2, whereas the azimuthal angle 
4>q varies over the course of the year, as shown in Fig. [2j If the earth's orbit were perfectly circular, then Q would 
increase at a constant rate, completing a full revolution in a year. In Sec. |III| we will see how to account for the fact 
that the earth's orbit is slightly elliptical, but for now we will simply express the azimuthal angle of the sun as a 
function of the time t. We therefore express the position of the sun, in the ecliptic frame of reference, as: 

/ cos </>©(£) \ 

r (i) = sin (t) . (2) 



B. Equatorial frame 

The axis of rotation of the earth is tilted with respect to the ecliptic frame by an angle of obliquity 

e = 23.44° = 0.4091 (3) 

(which corresponds to Milton's "twice ten degrees and more" PQ). It is therefore convenient to change coordinates 
to an "equatorial frame," by rotating about the rr-axis by an angle e, as shown in Fig. [3^a), so that the new z'-axis 
coincides with the earth's axis of rotation. The motion of the sun in this equatorial frame is illustrated in Fig.[3j[b), in 
which the celestial north pole is labelled P and the celestial south pole P. The ecliptic intersects the celestial equator 
at two points, e and e, known as the equinoxes. At e the sun crosses the equator from south to north, and this is 
therefore known as the northward equinox (or "first equinox", since it occurs first in the calendar year). Conversely, 
e is known as the southward, or second equinox. The points of maximum displacement between the position of the 
sun and the celestial equator are known as solstices, and are marked in Fig. plb) by s and s. 



2 The ecliptic is sometimes defined as the plane of the earth's orbit around the sun. The circle that we call the "ecliptic" is the intersection 
of that plane with the celestial sphere. 

3 For this reason the distant stars, which form the constellations, are also referred to as the "fixed stars." 
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FIG. 3: (a): We may transform from the ecliptic to the equatorial reference frame by rotating along the a;-axis by an angle 
equal to the obliquity e, given in Eq. so that the new axis z'-axis is also the axis of the earth's rotation, (b): The motion 
of the sun along the ecliptic, as seen in the new equatorial reference frame. The point P marks the celestial north pole and 
P the celestial south pole. The northward and southward equinoxes are marked by e and e, respectively. The northern and 
southern solstices are indicated by s and s, respectively. 



The position of sun in the equatorial coordinate frame is given by: 



The polar angle for the sun in this equatorial reference frame is therefore 

6q = arccos^Q = arccos (sine sin (/> ) . (5) 

In the astronomical literature, the value of measured with respect to the first equinox e, is called the sun's 
"ecliptic longitude." The corresponding 0q in the equatorial frame is called the sun's "right ascension," while tt/2 — 9'q 
is its "declination." For a full discussion of the celestial sphere and of the coordinate systems that astronomers use 
to characterize points on it, see [5]- 



sin 6q cos </>q 
sin O'q sin (f>' Q 
cos 8q 

1 

cos e — sin e 
sin e cos e 




cos <p 
cos e sin 0© 
sin e sin 4>q 



(4) 



C. Terrestrial frame 



Seen from the earth, objects in the sky rotate azimuthally in the equatorial frame (i.e., about the z'-axis), with 
constant angular velocity 



23.9345 hours ' v ; 

where 23.9345 hours is the duration of the "sidereal day," equal to the amount of time that it takes the earth to 
complete one rotation about its axis (and therefore also for a distant star to return to the same position in the sky). 
This is slightly less than the "mean solar day" of 24 hours, because of the sun's motion along the ecliptic (with respect 
to the distant stars) during the course of one sidereal day. 

To characterize the position of the sun, as seen from a point on the surface of the earth, we must also adjust for the 
geographic latitude L. We can achieve this by rotating about the x-axis by an angle equal to the co-latitude 7r/2 — L. 
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The transformation from the equatorial frame to the terrestrial frame therefore gives: 




1/0 

z Q ) \ cos<9 j 

cos — to)] sin [oj(t — t )] 0\ / cos</>© 
- sin [ui(t — to)] cos [oj(t — to)] cose sin ^ Q 
1 J \ sine sin 00 

where t — to is the interval during which the earth has rotated, measured with respect to a reference time 



(7) 



We 



will discuss how to choose the value of to (which will depend on the geographic longitude I) in Sec. VI 

By Eq. Q, the altitude (or "elevation") of the sun above the horizon, as a function of the latitude and the time t, 

is 



a & {L,t) = --6%{L,t) 



= arcsin[zQ(L,£)] 

= arcsin (— cos L ■ cos[</>0 (t)] • sin [uj(t — to)] + cos L ■ cos e • sin[0 (t)] ■ cos [to(t — to)] 

+ sin L ■ sine • sm[0©(t)]) . (8) 



When a = 0, the sun is either rising or setting. When a Q = 7r/2, the sun is directly overhead, at the "zenith" 
(which can occur only at tropical latitudes — e < L < e). 

Meanwhile, the azimuthal angle Q can be computed from Eq. (|7|), using the relation 



tan ( 



Ve 



(9) 



In Sec. VI we will work out the relation between this 4>q and the cardinal directions (North, East, South, and West). 



III. EQUATION OF THE CENTER 



For some purposes, it may be acceptable to approximate the angle <f>@{t) as increasing linearly and completing a 
full revolution in one year. A more precise expression can be obtained from Kepler's first and second laws of planetary 
motion, which state that the earth moves along an ellipse, with the sun at a focus, while the line segment from the 
sun to the earth sweeps out equal areas in equal times. 

The angle subtended by the line from the sun to the earth, with respect to the major axis of the elliptical orbit, 
is known to astronomers as the "true anomaly" and is usually represented by the letter v. Finding v as a function 
of time has no exact analytic solution, 5 but an expansion can be obtained, which converges rapidly for small orbital 
eccentricity e, known as the "equation of the center:" 

w = M + 2esinAf + -e 2 sin2M+ — e 3 (13sin3M - 3sinM) + C(e 4 ) . (10) 
4 12 



The "mean anomaly" in Eq. ( 10 1 can be expressed as 

M = Mo + M x t , (11) 
with constant Mo,i; it would be equal to the angle v for a perfectly circular orbit (e = 0) of equal area to the true 



elliptical orbit; see [8]. The values of v and M in Eq. (10) are measured with respect to the perihelion, which is the 
point of closest approach between the earth and the sun, as shown in Fig. |4j The value of 2-k/M\ is slightly greater 
than one calendar year because of the slow precessions of the equinoxes and the perihelion, which we will discuss in 
Sec. HITaI 



4 Note that the signs of the ± sin[o!(t — to)] off-diagonal entries in the corresponding rotation matrix in Eq. |7| reflect the fact that the 
earth's rotation displaces the sun in an azimuthal direction opposite to that of the sun's yearly motion along the ecliptic. This is the 
reason why the mean solar day of 24 hours is longer than the sidereal day of 23.9345 hours: the extra 4 minutes of rotation are need to 
compensate for the change in 00 in Eq. (|2j. 

5 Newton offered a rigorous proof that no such analytic solution could exist, using concepts now associated with topology, long before 
topology was invented. This fascinating proof (the first impossibility proof since the ancient Greeks) is discussed in [6l |7|. 
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FIG. 4: Diagram of the earth's orbit around the sun. The eccentricity is exaggerated for clarity. The "true anomaly" v 
is measured from perihelion, which is the point of closest approach between the sun and the earth. (The point of greatest 
separation between sun and earth is called the "aphelion.") We measure the ecliptic azimuthal angle of the sun, <f)Q, from the 
position at the first equinox. Therefore (j>® = v — Vq, where vo is the angular displacement between the perihelion and the 
equinox. 

For our purposes it will be convenient to measure the angle Q from the first equinox of the year. Therefore we let 

00 = « - v , (12) 

where vq is the angular displacement between the perihelion and the first equinox, as shown in Fig. [4] Using the 
current astronomical data for the parameters Mo, Mi, e, and vq [5], we can write the equation of the center for the 
earth as: 

M(t) = -0.0410 + 0.017202 1 (13) 

and 

(t) = -1.3411 + M(t) + 0.0334 sin[M(t)] + 0.0003 sin[2M(t)] , (14) 

where t = corresponds to 1 January 2013, 0:00, Universal Coordinated Time (UTC), and t is measured in mean 
solar days of 24 hours. 

A. Precession of equinoxes and perihelion 

In the second century BCE, the Greek astronomer Hipparchos of Nicaea found that the positions of the equinoxes 
moved along the ecliptic (i.e., with respect to the distant stars) by about 1° per century (the modern estimate is 1.38° 
per century) . Newton correctly explained this as due to the tidal forces that the moon and the sun exert on the earth, 
which is not perfectly spherical. If the earth did not spin, those tidal forces would pull the earth's equatorial bulge 
onto the orbital plane of the corresponding perturbing body (i.e., of either the moon or the sun). The earth's spinning 
turns the action of that tidal torque into a precession, so that the axis of the earth's rotation describes a cone, and 
the position of the celestial north pole therefore moves slowly along a circle, with respect to the constellations. 6 



6 This slow change of the positions of the poles, equinoxes, and solstices, relative to the distant stars, implies that the signs of the Zodiac 
are not fixed with respect to the solar calendar. For example, the "Tropic of Cancer" (Milton's "Tropic Crab" in [T]) was so named 
because the position of the sun at the time of the northern solstice used to lie within the constellation of Cancer, but today the northern 
solstice actually lies in Taurus. The first equinox, which used to lie in Aries when the ancient Babylonians developed the calendar, has 
since shifted to Pisces and will move into Aquarius around the year 2,600. This last circumstance has been the source of much mystical 
twaddle about the "dawning of the Age of Aquarius." 




FIG. 5: (a): Cross-section of the celestial sphere along the earth's axis of rotation PP, centered at the position O of an 
observer at geographic latitude L. The point m corresponds to the maximum altitude of the sun, and m to the minimum 
altitude, (b): Cross-section of the celestial sphere, centered at point a and perpendicular to the earth's axis of rotation. The 
point c corresponds to sunrise and c to sunset. The arrows show the direction in which the celestial sphere rotates with respect 
to the observer at O. 



The period of the precession of the earth's axis is about 26,000 years. Since the recurrence of the seasons depends 
on the periodicity of the equinoxes, rather than on the actual time it takes the earth to go once around the sun, the 
modern calendar is based on the "mean tropical year," which is shorter than the sidereal year by about 20 minutes 
(i.e., 1/26,000 of a sidereal year). 

The position of the perihelion with respect to the distant stars also varies, but more slowly, with a period of 
about 112,000 years, which is equivalent to a displacement of about 0.32° per century. This precession results from 
perturbations to the motion of the earth around the sun caused by the gravitational pull of the moon and the other 
planets, and to a lesser extent also by relativistic corrections to Newtonian gravity. 7 

The respective precessions o f th e equinox and that of perihelion proceed in opposite directions along the ecliptic, 
causing the value of Vq in Eq. (12) to decrease by about 1.7° per century. 8 Though for our purposes such precision is 
hardly justified, if we wished to take into account the precession of the equinoxes and the perihelion, we could make 
vq in Eq. ( 12 ) a time-dependent parameter. 



IV. DURATION OF DAYLIGHT 



The computation only up to Eq. ^ suffices to obtain a good estimate of the number of hours of daylight for a given 
day of the year, if we do not care for the precise time of sunrise and sunset. Here the main approximation is that 
that the azimuthal angle of the sun in the ecliptic frame, 4>q, will be taken to be fixed during a given calendar date 
d. For definiteness, let us say that 4>q is computed at noon for the date and location of interest, the corresponding 
time being translated to Universal Coordinated Time (UTC) for use in Eqs. (14) and (13). 

As the sphere 



Figure [5[a) shows a cross-section of the celestial sphere, parallel to the earth's axis of rotation PP. 
rotates about the observer at point O, the celestial pole P maintains a fixed altitude, equal to the observer's geographic 



7 One of the most convincing early demonstrations of the validity of Einstein's theory of general relativity was that it explained the 
anomalous precession of the perihelion for the orbit of Mercury, which astronomers had until then failed to account for by the gravitational 
influence of the known planets. This subject is treated in detail in |10j . 

8 The quantity 2tt — vq is known to astronomers the "longitude of perihelion." 
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FIG. 6: N umb er of hours of continuous daylight H, as a function of the day of the year d (starting on 1 January), computed 
using Eq. ( |17| , for: the latitude of Cartagena de Indias, Colombia, 10°24' N (red curve); the latitude of Boston, Massachusetts, 
USA, 42°21 7 N (blue curve); the latitude of Stockholm, Sweden, 59°20' N (green curve); and the Arctic Circle, 66°34' N (dashed 
black curve). 



latitude L. 9 Point m marks the maximum altitude of the sun, while point m marks its minimum altitude. 

The path of the sun in the sky corresponds to the circle am, shown in Fig. [5jb) (again, as long as we neglect the 
change in <f>Q, and therefore also in 9'q, during the course of one day). This circle is a cross-section of the celestial 
sphere, perpendicular to the axis PP and parallel to the line mm. 

In terms of the angle 8 in Fig. p[b), 10 the number of hours of daylight is simply 



if = 24 1 



(15) 



since the sun moves along the circle am uniformly, with a period of 24 hours. 11 Examining Figs.|5|a) and (b), we see 
that 



ab 



5 = arccos = arccos (tan L cot 9' Q ) 



am 



(16) 



Therefore we can express the number of hours of daylight as a function of geographic latitude and day of the year in 
the form: 



H(L,d) = 24 



= 24 



arccos (tan L cot 9' Q (d)) 

7T 



1 arccos | tan L 

7T 



sine sin < 



v/l — sin 2 e sin 2 [4> Q (d)] t 



(17) 



Figure [6] shows plots of if as a function of the day of the year d, at the latitudes of Cartagena de Indias (Colombia), 
Boston (USA), Stockholm (Sweden), and the Arctic Circle, all in the northern hemisphere. Note that, for the northern 



9 If we take L to be positive for points on the northern hemisphere of the earth, then P is the north celestial pole, and P is the south 
celestial pole. The opposite convention would be more convenient for observers in the southern hemisphere. 

An astronomer would call <5 the sun's "local hour angle" at the times of rising and setting. See, e.g., 

1 By making the period of rotation of the sun about the celestial poles in Fig.|5|equal to the mean solar day of 24 hours, rather than the 
sidereal day of 23.9345 hours, we are taking into account the average azimuthal displacement of the sun during the course of one day. 
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FIG. 7: The upper red curve corresponds to the maximum altitude of the sun, a m , in degrees, as a function of the day of 
the year d (starting on 1 January), at the latitude of San Jose, Costa Rica (9° 56' N). The lower red curve gives the minimum 
altitude a™ at that same latitude. The dashed blue curves give a m and a™ for the latitude of Casablanca, Morocco (33° 32' 
N). 



hemisphere, the midyear solstice (which occurs around 21 June, or d = 171) is always the longest day, whereas it is 
the shortest day everywhere in the southern hemisphere. Conversely, the year-end solstice (around 21 December, or 
d = 354) is always the longest day in the southern hemisphere and the shortest in the northern hemisphere. 



A. Maximum and minimum solar altitudes 



In Fig. [5] it is easy to see that the maximum and minimum altitudes of the sun on a given date, which we respectively 
are: 



label a rn and a m , 



a m (L,d) — arcsin (sin \L + 6' & (d)~\ ) 
am(L,d) — arcsin (sin [L — 0q(c?)]) , 



(18) 



where 6L is given by Eq. (l5|) . These are the values between which the solar altitude of Eq. (l8|) varies during the day 



The use of the arcsine function in Eq. (18) reflects the fact that an altitude must, by definition, be between ±7r/2. 

For tropical latitudes — e < L < e it is not true that the sun reaches its maximum altitude over the horizon on the 
day of the solstice. This is shown graphically in Fig. [7| which plots the maximum and minimum solar altitudes, given 
by Eq. (18), as functions of the day of the year, for the latitudes of San Jose, Costa Rica (9°56' N) and Casablanca, 
Morocco[33 32' N). While the noonday sun over Casablanca does reach its maximum altitude (ir/2 — L + e) on the 
day of the midyear solstice, the sun over tropical San Jose reaches zenith (a m = n/2) on two days, one before and 
the other after the solstice. We shall compute these dates in Sec. VII E 



B. Adjustments 



Comparison of the results of Eq. (17) with the hours of sunrise and sunset published in newspapers and other 
sources reveals a discrepancy: we underestimate the duration of daylight by several minutes. One reason for this 
discrepancy is that the times of sunrise and sunset are usually taken to correspond to the moments when the upper 
edge of the solar disk crosses the horizon. Since the solar disk has an angular diameter of about 0.5°, sunrise occurs 
slightly before, and sunset slightly after, the times that we have computed, which referred to the horizon crossings of 
the center of the sun. 

An even more important factor is that, since the density of the earth's atmosphere increases closer to the earth's 
surface, light passing obliquely through the atmosphere bends downwards. This atmospheric refraction appreciably 



11 



c'/ 

/ 1 
/ 1 






/ 1 
I 1 


b 




1 
1 

\ 1 


j 


V 


\ 1 
c'\ 







(a) 




FIG. 8: (a): Cross-section of the celestial sphere, centered at point a and perpendicular to the earth's axis of rotation PP; see 
Fig. [SJb). The points c' and c' correspond to sunrise and sunset, respectively, after taking into account the size of the solar 
disk and the effect of atmospheric refraction, (b): Cross-section of the celestial sphere, along the earth's axis of rotation and 
centered at the observer's position O; see Fig. pTa). The vertical displacement between b and b' corresponds to the terrestrial 
coordinate of the sun z'q (see Eq. Q) at the time of sunrise or sunset. 



lifts the sun's apparent position when it is near the horizon. When the sun (or any other celestial object) appears to 
us to be at horizon, its true altitude is about —0.6° [TTj . 

In [12] , Saemundsson proposed a simple formula for the apparent shift in altitude, Aa, as a function of the true 
altitude a, accurate to within 0.07' for all values of a: 



Aa 



1.02' 



tan a 



10.3° 
a+5.11 c 



X 1° 



(19) 



This expression may be used to adjust the altitude of Eq. Q in order to bring it into closer agreement with observation. 
(For a full treatment of atmospheric refraction and other effects on the apparent position of celestial objects, as seen 
from the surface of the earth, see [15] •) 

Combining the size of the solar disk with the lifting due to refraction, we see that sunrise and sunset correspond 
to the moments when the true altitude of the center of the solar disk is a© = — (0.5°/2 + 0.6°) = —0.85°. Thus, it is 
necessary to adjust the diagrams of Fig. [5] to account for this displacement, in order to bring our expression for the 
duration of daylight into agreement with the published times for sunrise and sunset. 

Points c' and c' in Fig. [8]ja) mark the position of the sun at the published times of sunrise and sunset, respectively. 
Point b 1 marks the projection of these unto the line mm. The correction to the duration of daylight, in hours, is 
therefore 



AH = 24 



5-5' 



(20) 



In order to obtain an expression for the adjustment AH we therefore need to express 5 — 5' as a function of the 
latitude and the date of the year. 

The vertical direction in Fig. Isfb) corresponds to z"-axis in the terrestrial reference frame of Eq. ([7]). The vertical 
displacement between b and b' therefore corresponds to the value of z'q at the published times of sunrise or sunset, so 
that 



sin 0.85° 



0.015 



cos L cos L cos L 
Meanwhile, we can see from the diagrams in Fig. [8] that 

bb' = ab' — ab = sin Q' Q (cos 5' — cos 5) ~ sin 6L sin 5(5 — 5') 



(21) 



(22) 
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FIG. 9: Cross-section of the celestial sphere — as shown in Fig. [5|a) — along the horizontal plane passing through the position 
of the observer at O. The ray from O to b points North for an observer in the northern hemisphere (South for an observer in 
the southern hemisphere). The sun rises at c and sets at c. The angle ip c therefore corresponds to the geographic azimuth (i.e., 
compass bearing) of sunrise. 



Combining Eqs. (20), (21), and (22), we obtain 



AH 



7 min. 



cos L sin 6'q sin S 



(23) 



We can express the value of AH as a function of latitude L and date d by plugging in the expression for 9' & in Eq. (p>|) 
and for S in Eq. (16). Equation (23) implies that the adjustment to the duration of daylight from the size of the 



solar disk and atmospheric refraction amounts to at least 7 minutes, for all locations and dates of the year, and can 
be significantly larger for locations distant from the Equator (L — 0) and on dates far from the equinoxes (when 
6> q = 6 = k/2)P 

Note also that the position of the sun with respect to the horizon is further lifted, and the duration of daylight 
consequently lengthened, by observing from an elevated position, since this makes the viewer's horizon recede. Indeed, 
it is possible to estimate the size of the earth from the time difference between the sunset over the ocean as seen by 
an observer lying on the ground and the sunset as seen by the same observer standing up; see [14] . 



V. DIRECTION TO SUNRISE AND SUNSET 



Knowing 9L, we can also compute the direction, with respect to the cardinal points, in which an observer will see 
the sun rise and set. Figure [9] shows a cross-section of the celestial sphere, this time through the horizontal plane 
passing through the position of the observer at O, so that the resulting circle has the same unit radius as the celestial 
sphere itself. If the observer is on the northern hemisphere, then the ray from O to b points North (whereas it points 
South for an observer in the southern hemisphere). 

The angle ip c therefore gives the geographic azimuth (i.e., the compass bearing) of sunrise and can be expressed as: 

/ cos 8'p 

w c = arccos Ob = arccos - 

\ cosL 

If we neglect the change in the position of the sun between sunrise and sunset on the same day, then the azimuth of 
the sunset at c is simply ip s = 2-7T — tp c . 




Equation |23| should not be used at Arctic or Antarctic latitudes (i.e., \L\ > it/2 — e), or very near them, because (<5 — S')/S might not 
always be small, invalidating the approximation in Eq. {22} , and also because o m and may have the same sign, implying that the 
sun neither rises nor sets. 
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FIG. 10: Plot of the geographic azimuth for sunset, (pc, in degrees, at the latitude of Manhattan, in New York City (40°47' 
N), as a function of the day of the year d (starting on 1 January). When (p B = 299° (marked by the dashed line) the sunset is 
aligned with the east-west streets on the main traffic grid. 



A. Manhattanhenge 



Figure 10 gives a plot of the sunset's azimuth tp E , as a function of the date d, for the latitude of Manhattan. This 
can be used to find the dates of "Manhattanhenge" (also called the "Manhattan solstice"), when a New York City 
pedestrian can see the sunset in between the skyscrapers, because the sunset is aligned with the east-west streets on 





FIG. 11: Photograph taken by the author on 11 July 2012, at 8:21 p.m., Eastern Standard Time (EST), looking west along 
34th Street, from a position near the intersection with Park Avenue, in midtown Manhattan, New York City. The solar disk 
can be seen between the profiles of the buildings, a few degrees above the horizon. The Empire State Building is prominent on 
the left. 
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the main traffic grid for the borough of Manhattan [13] ■ Since those streets point 29° north of the true west, this 
alignment occurs when ip E (d) = 270° + 29° = 299°. For the year 2013 these correspond to 28 May (d = 147) and 14 
July (d= 194). 

A more impressive visual spectacle than the alignment of the actual sunset with the east-west streets is the ob- 
servation of the full solar disk in between the profiles of the buildings and slightly above the horizon. This occurs a 
couple of days after the first date of the year for which (p E (d) — 299°, as well as a couple of days before the second 



such date (see [H]). Figure 11 shows a photograph of the view on 11 July 2012, looking west along 34th Street in 
midtown Manhattan, several minutes before sunset. 

The name Manhattanhenge refers to the fact that Stonehenge, in the English county of Wiltshire, was built in 
prehistoric times in such a way that around the time of the summer solstice an observer standing at the center of the 
circular structure sees the sun rise above the outer Heelstone. There are other famous instances of ancient structures 
oriented according to solar alignments. The Great Temple of Amen-Ra, at Karnak, Egypt, was designed so that 
the last rays of the sun on the day of the summer solstice illuminate the inner sanctuary. 13 The Mayan pyramid of 
Kukulkan (also called El Castillo, "The Castle" ) in Chichen Itza, in the Mexican state of Yucatan, was built in such 
a way that the setting sun, around the time of the equinoxes, casts a shadow that looks like a serpent slithering down 
the side of the staircase. 



VI. ACCOUNTING FOR GEOGRAPHIC LONGITUDE 



There is one issue left to resolve before we can give the sun's position in the sky for a given time and location: the 
choice of to in Eq. (l8|) , needed to translate from the 4>q in Eq. (l7|) to a geographic azimuth ip & defined with respect to 
the cardinal points (North, East, South, and West). This must De consistent with how t is measured in Eq. (13). The 



choice of to will evidently depend on the observer's geographic longitude, since at a given universal time the azimuth 
of the sun with respect to the cardinal points depends on the longitude at that location. 

At the time of the northward equinox, when 4>q — 0, the sun's altitude, observed at the Equator (L = 0), is given 
by Eq. ^ as a = lo • (to — t). Therefore, we can simply find the geographic longitude lo at which the equinox occurs 
precisely at sunset and, for any other longitude £, choose: 

W) = icq - '—^ , (25) 
LO 

where t eq is the time of the northward equinox and lo is the rate of the earth's rotation (see Eq. ([6])). The longitude 
I is taken to be positive to the east of the Prime Meridian. 

For instance, the first equinox for the year 2013 occurs on 20 March at 11:02 UTC [TBJ, so that t eq = 78.46 days. At 
that moment, the difference between apparent solar time and the "mean solar time" (which matches UTC for I = 0) 
is about 7 minutes |17j . Thus, the mean solar time for sunset at the time of the equinox is 18:07. 1 The first equinox 
of 2013 therefore occurs at sunset at the geographic longitude of 

i Q.n7 ii -no 

to = 2Am x 2 ^ - L854 = 106 - 3 ° > ( 26 ) 

to the east of the Prime Meridian. 



A. Geographic azimuth 



Note that </>q decreases during the course of a day, and that we have chosen to in Eq. ([7]) so that <f>o(L, lo, t cq ) = 0. 
The sunset at the time of the equinox must point directly West. Conventionally, the geographic azimuth tp is defined 
so that ip = corresponds to North, tp — 90° to East, ip = 180° to South, and ip = 270° to West. Therefore we should 
take 

p & = 270° - ^ . (27) 



This might remind some readers of the method by which Dr. Indiana Jones locates the Well of Souls in the movie Raiders of the Lost 
Ark (1981). 

The variation in the difference between apparent solar time and mean solar time throughout the year is given by the "equation of time." 
We shall return to this issue in Sec. IVIIII 



15 



VII. ALTITUDES 



We shall now plot the altitude of the sun above the horizon as a function of time (measured over the course of one 
year), for several geographic locations of interest. These are the results obtained from Eqs. (14), (13), and (25). 
for the center of the solar disk. 15 



A. Buenos Aires, Argentina 

Buenos Aires, the capital of Argentina, is located a latitude 34°36' South, longitude 58°23' West. The pattern of 
the sun's altitude here is typical of southern temperate regions, with longer days around the beginning and end of the 
calendar year, and shorter days around midyear. 
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FIG. 12: Altitude of the sun above the horizon, qq, in degrees, at the coordinates of Buenos Aires, Argentina (34°36' S, 58°23' 
W), as a function of the time of the year t (measured in days from 1 January, 2013, 0:00 UTC). 



B. Reykjavik, Iceland 

Reykjavik, the capital of Iceland, is located at latitude 64°08' North, longitude 21°56' West. Being south of the 
Arctic Circle (i.e., latitude 90° — e = 66°34'), the sun rises and sets over Reykjavik every day of year, but around the 
summer solstice, in late June, it never dips far enough below the horizon for the sky to grow completely dark. 
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FIG. 13: Altitude of the sun above the horizon, qq, in degrees, at the coordinates of Reykjavik, Iceland (64°08' N, 21°56' W), 
as a function of the time of the year t (measured in days from 1 January, 2013, 0:00 UTC). 



15 These plots were made without using the correction for atmospheric refraction of Eq. {l9| , which in any case is far too small an 
adjustment to be discernible at the resolution of our graphs. 
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C. Alert, Nunavut, Canada 



Alert, in the Canadian territory of Nunavut, is the northernmost permanently inhabited place on earth, located 
817 km from the geographic north pole. Its geographic coordinates are latitude 82°30' North, longitude 62°20' West. 
Here the sun rises and sets only in March and April, and again in September and October. During the rest of the 
year the sun remains either above or below the horizon. 




FIG. 14: Altitude of the sun above the horizon, ag, in degrees, at the coordinates of Alert, Nunavut, Canada (82°30' N, 62°20' 
W), as a function of the time of the year t (measured in days from 1 January, 2013, 0:00 UTC). 



D. Singapore 

Singapore, a city that is coextensive with the independent nation of the same name, is located very near the 
Equator, at latitude 1°17' North, longitude 103°50' East. Here the days remain almost evenly divided between light 
and darkness. The sun only reaches zenith at around the times of the equinoxes (in late June and late September). 




FIG. 15: Altitude of the sun above the horizon, a©, in degrees, at the coordinates of Singapore (1°17' N, 103°50' E), as a 
function of the time of the year t (measured in days from 1 January, 2013, 0:00 UTC). 



E. San Jose, Costa Rica 



San Jose, the capital of Costa Rica, is located at latitude 9° 56' North, longitude 84° 5' West. The pattern of the 
sun's altitude is typical of northern tropical regions, with slightly longer days around the midyear solstice. Note that 
the noonday sun reaches the zenith around two different dates, which by Eq. ( 18 1 are given by 



9'(d) = 90° - 9°56' = 80.07° = 1.397 . 



(28) 




FIG. 16: Altitude of the sun above the horizon, qq, in degrees, at the coordinates of San Jose, Costa Rica (9°56' N, 84°5' W), 
as a function of the time of the year t (measured in days from 1 January, 2013, 0:00 UTC). 



VIII. ANALEMMAS AND EQUATION OF TIME 



The mean solar time is defined in terms of a fictitious "mean sun," which moves at a uniform rate along the celestial 
equator. This is different from the actual position of the sun, which moves along the ecliptic at a rate which is not 
quite uniform (see, e.g., |18|). At noon, the fictitious mean sun crosses the meridian above the observer. Here, the 
word "meridian" is understood not as a line of constant geographic longitude, but rather as the great circle OP in 
Fig. [5^ a), which is perpendicular to the observer's horizon and passes through the celestial poles. The hour of 12:00 
UTC corresponds to the time when the mean sun, as seen by an observer at any location along the Prime Meridian 
(longitude £ = 0), has a positive altitude and a geographic azimuth of either <p = 180° (for observers in the northern 
hemisphere) or ip — (for observers in southern hemisphere). 

If we plot the actual sun's altitude a© vs. its azimuth ip@, over the course of many days, for a given mean solar 
time and at a given location, the result will look like a figure 8, called an analemma. Figure [17| illustrates this for the 
geographic coordinates of Greenwich, England (51°29' N, 0° E). These plots correspond, respectively, to fixing the 
hour of observation at 0:00, 6:00, 12:00, and 18:00 UTC, for each day of the year 2013. 

The difference, as a function of the date of the year, between noon and the actual time when the sun crosses the 
meridian, is called the "equation of time." This difference results mainly from the projection of points along the 
ecliptic onto the equatorial plane (i.e., from the coordinate transformation between the ecliptic and equatorial frames, 
given by Eq. (HI), which implies that tan^Q = cose tan </>© and therefore distorts 0q with respect to 0©), and to a 
lesser extent also from the eccentricity of the earth's orb it, which makes the motion of the sun along the ecliptic, 
described by 0©, less than perfectly uniform (see Sec. Ill I. 16 

The equation of time can be read from the horizontal displacement of the analemma in Fig. |17[ c). In that plot, 
each degree of azimuthal displacement away from (p® = 180° translates to 24 hours /360 = 4 minutes in the equation 
of time for the corresponding date. 

Even though the position of the sun in the sky for a given hour, calendar date, and location, is not quite constant 
from one year to the next, the way in which the mean solar time is defined ensures that the overall shape of the 
analemma remains unchanged. The word "analemma" derives from the Greek term for the pedestal of a sundial and 
refers to the instructions on how to correct a sundial's reading for the equation of time, which can give an adjustment 



Introductory textbooks and other pedagogical sources sometimes explain the analemma's horizontal displacement as resulting solely or 
primarily from the eccentricity of the earth's orbit (see, e.g., I19p . This is clearly incorrect: Kepler's second law implies that a planet's 
orbital angular velocity is faster than the mean when the planet is near the perihelion and slower when it is near the aphelion (see 
Fig. [4). If this were the dominant factor in determining the equation of time, then the azimuth of the sun, observed at noon, would 
be greater than 180° during half of the year and smaller than 180° during the other half, giving a figure rather than a figure 8 for 
the analemma. Though the orbital eccentricity does play a significant role in determining the precise form of the equation of time, the 
analemma would still be a figure 8 if the earth's orbit were perfectly circular. 
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FIG. 17: Plots of solar altitude oq vs. azimuth 950, in degrees, for each day of 2013, at the geographic coordinates of Greenwich, 
England (51°29' N, 0° E), for the hours of: (a) 0:00, (b) 6:00, (c) 12:00, and (d) 18:00 UTC. 



of as much as a quarter of an hour. For a full treatment of these issues as they relate to the construction and reading 
of sundials, see [20] . 
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